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Abstract. We propose a simple conjecture for the functional form of the asymptotic behavior of work 
distributions for driven overdamped Brownian motion of a particle in confining potentials. This conjecture 
is motivated by the fact that these functional forms are independent of the velocity of the driving for all 
potentials and protocols, where explicit analytical solutions for the work distributions have been derived in 
the literature. To test the conjecture, we use Brownian dynamics simulations and a recent theory developed 
by Engel and Nickelsen (EN theory), which is based on the contraction principle of large deviation theory. 

Our tests suggest that the conjecture is valid for potentials with a confinement equal to or weaker than 
the parabolic one, both for equilibrium and for nonequilibrium distributions of the initial particle position. 

For potentials with stronger confinement, the conjecture fails and gives a good approximate description 
only for fast driving. In addition we obtain a new analytical solution for the asymptotic behavior of 
the work distribution for the V-potential by application of the EN theory, and we extend this theory to 
nonequilibrated initial particle positions. 

PACS. 05.40.-a Fluctuation phenomena, random processes, noise, and Brownian motion - 05.70.Ln 
Nonequilibrium and irreversible thermodynamics 

1 Introduction 

With the ever-improving possibilities to manipulate and 
control physical systems on the single-molecule level, the 
analysis of the role of thermodynamic quantities such as 
work, heat, and entropy, when defined for single system 
trajectories, has become an important field in nonequilib¬ 
rium statistical mechanics mm- One of the most intrigu¬ 
ing achievements in this field is the discovery of detailed 
and integral fluctuation theorems [3J21E1EHZ1E], which in 
general refer to nonequilibrium systems that are externally 
driven by a protocol A (t) of control variables during a time 
interval [0 ,tf]. These theorems hold true universally and 
they can be viewed as a generalization of the second law 
of thermodynamics. 

Among the theorems, the detailed Crooks fluctuation 
theorem (CFT) [T], and its integral counterpart, the Jarzyn- 
ski equality (JE) [3] are perhaps the most prominent. 

They pertain to systems in contact with a heat reservoir 
at temperature T and in equilibrium at the initial time 
t a = 0. The CFT states that work probability distribu¬ 
tions (WPD) p(w) and pr(w) for a protocol A(i) and the 
associated reversed protocol A n(t) = A (tf — t), respec¬ 
tively, are related according to p(w)/pr(—w) = eP( w ~ AF ") 
g], where AF = Ff — Fq is the difference in equilib¬ 
rium free energies Fqj of the macrostates specified by 


the control variables Xqj = X(t 0 j), and f3 = (k-eT)^ 1 is 
the inverse thermal energy (fep is the Boltzmann constant 
and T is the temperature). The JE states that (e~^ w ) = 
f dw p(w)e~ /3w = e~ /3AF [3|. It becomes particularly valu¬ 
able in unidirectional experimental settings, where the 
work for the reversed protocol can not be measured, and 
accordingly the Crooks theorem can not be applied [5]. 
In applications of the JE, measured histograms generally 
need to be extended to the tail regime, because the aver¬ 
age of e~ /3w in the JE is dominated by rare trajectories 
with work values w <C AF. This problem can in principle 
be resolved by fitting the wings of a measured WPD to 
theoretical predictions for the asymptotic behavior in the 
limit w —> — oo. 

Corresponding predictions for this asymptotic behav¬ 
ior are, however, difficult to obtain. WPD depend on de¬ 
tails of the experimental setup and only a few generic 
properties have been reported so far. In the limit of quasi¬ 
static driving, the WPD becomes a delta function, p(W) = 
S(W — AF ), and for sufficiently slow driving, close to 
the quasi-static limit, it can be approximated by a Gaus¬ 
sian distribution in the relevant regime of the integrand 
in the JE m- The JE equality moreover gives the con¬ 
straint that p(w)e~ l3w is integrable, which implies that for 
w —> — oo, the WPD must decay faster than 
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For overdamped Brownian motion of a particle in time- 
varying potentials, a sophisticated theory based on the 
contraction principle of large deviation theory was recently 
developed by Engel and Nickelsen (EN theory) |11II12| . 
This allows one to predict the asymptotic behavior of the 
WPD by solving a system of ordinary differential equa¬ 
tions with certain boundary conditions (see Appendix [All. 
Analytical solutions of these equations were given for the 
harmonic potential mm, where in the protocol either 
its minimum is moved (“sliding parabola”) or the stiff¬ 
ness is varied (“breathing parabola”). In the Appendix iBl 
we further derive an analytical solution for the case of a 
V-potential |13] with time-varying slope. However, in gen¬ 
eral, the explicit functional form of the WPD asymptotic 
behavior remains unknown and the respective differential 
equations have to be solved numerically. 

In this work we show that the functional form can be 
often guessed by a simple method, which becomes exact 
in the limit of infinitely fast driving. It is referred to as 
the functional form (FF) conjecture in the following. We 
show in Sec. [2] that for the rare cases of potentials and 
protocols, where the WPD has been derived analytically, 
the FF conjecture always provides the correct functional 
form. This leads to the question whether the FF conjecture 
can be applied in general to predict the functional form of 
the WPD asymptotics. 

To tackle this problem, one could test the FF conjec¬ 
ture against simulations. However, in corresponding sim¬ 
ulations it becomes very difficult to capture the WPD be¬ 
havior for large negative work values with sufficient statis¬ 
tics. Another possibility is to compare the predictions of 
the FF conjecture with results obtained from the EN the¬ 
ory. This theory, which relies on the applicability of large 
deviation theory to the WPD, has been shown to repro¬ 
duce the exact Gaussian WPD for the sliding parabola El 
m and further results obtained from it for the breathing 
parabola turned out to agree with an exact treatment E- 
Moreover, the EN theory was shown to match simulated 
data for the asymptotic WPD behavior for several other 
potentials and protocols EE- 

Using the exact results for the logarithmic-harmonic 
potential m, we provide in Sec. [2] further severe evidence 
of the validity of the EN theory. Given this evidence, we 
then base our evaluation of the FF conjecture on the EN 
theory. This evaluation is supplemented by WPD data 
from Brownian dynamics (BD) simulations. In the corre¬ 
sponding analysis of the FF conjecture in Sec. [3l we will 
consider both “hard potentials” (with stronger confine¬ 
ment than the harmonic one) and “soft potentials” (with 
weaker confinement than the harmonic one). In Sec. [I] we 
furthermore generalize both the EN theory and the FF 
conjecture to the case of a non-equilibrated initial posi¬ 
tion of the particle. 


Table 1 . Tail behavior predicted by the FF conjecture, Eq. ©, 
for the potentials, where the (whole) WPD has been calculated 
analytically [for certain protocols A (t), e.g. the form in Eq. ©], 
and for the V-potential (last line), where the asymptotic be¬ 
havior of the WPD can be derived exactly (see Appendix [Bj). 
The parameters k > 0 and g > — /3 _1 in the potential forms 
in the first and third line are constants. The parameters A, B 
and C entering the functional form of the WPD tails are con¬ 
stants that depend on details of the protocol and the inverse 
temperature /3. 


Potential U(x, X(t)) 

Tail behavior 

f[*-A W] 2 

Ae-^-c) 2 

|A(t)a; 2 

A\w\- 1/2 e~ B ^ 

~9 log 1*1 + 5 A (t)x 2 

/l|-u;| _(1 -'" )/2 e Hw 

X(t)\x\ 

Ae“ sw 


potential is given by 


dx dU(x,X(t)) 
dt ^ dx 


+ pit )> 


( 1 ) 


where p is the mobility, U{x,X(t)) is the potential, and 
rj[t) is a Gaussian white noise with zero mean and corre¬ 
lation {p{t)p{t')) = 2pkB,T5{t — t'). In the following we set 
k^T as our energy unit, and consider x as dimensionless, 
which allows to set {pk^T) -1 as the time unit [16] . The 
work w done on the particle along the stochastic trajec¬ 
tory x(t) during the time interval [0,f/] is 


w = 


f tf dt dU(x(t),X(t))d\(t) 

Jo 9X dt 


( 2 ) 


To keep the treatment simple, we consider potentials that 
depend monotonically on the protocol A(f) (the only ex¬ 
ception is the sliding parabola potential in the first line 
in Tab. [[]), and we always use monotonic protocols, where 
both dU(x(t),X{t))/dX and dX{t)/dt do not change sign 
in [0,t f ]. 

In the limit of infinitely fast driving, where the protocol 
A(f) jumps from its initial value Ao = A(0) to its final value 
Xf = X(tf) at some time instant in the interval [0,f/], the 
WPD is given by p(w) = f dx S [w — AU{x)] po(x) where 
AU{x) = U{x , A (tf)) — U(x , A(0)), and po{x) is the initial 
distribution of the position of the particle at time t = 0. 
In Refs. pfT7| it was pointed out that for the breathing 
parabola the functional form of the WPD asymptotics is 
the same as that for infinitely fast driving. 

This suggests the following FF conjecture: If the po¬ 
sition of the Brownian particle is initially in equilibrium 
with distribution p eq (x) = exp [— (3U(x, A(0))]/Zo, the func¬ 
tional form of the WPD asymptotic behavior is given by 


2 FF conjecture and EN theory versus exact 
results 


p(w) = 

tun 


dx 6 [w ~ AU(x)\ Pe q(x) , w— >-± 00 . 


(3) 


Here the symbol “ = ” indicates that the left and the right 

J fun 6 


The time evolution of the particle position x for over- hand-side equal in their functional form, but may dif- 
damped one-dimensional Brownian motion in a time-varying fer in specific values of parameters. This means that in 
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Fig. 1. Asymptotic behavior of the WPD for the logarithmic-harmonic potential calculated from the EN theory (symbols) 
and fitted to the asymptotic behavior p(w) ~ |u>| - " exp(—£>|ia|) (lines) for protocol parameters r = tf = k = 1 [see Eq. (|4|)[. 
and different strengths (a) g = 1.5, (b) g = 10, and (c) g = 20 of the logarithmic part. The insets show the convergence of 
the fitting parameters v and B to their correct values B ~ 2.335 (independent of g) and v = (1 — (3g)/2 (cf. third line in 
Tab. [TJ with decreasing upper bound of the fitting interval [—320, ui+]. The errors bars given in the insets, which are smaller 
than or comparable to the size of the symbols, mark 95% confidence bounds of fitting parameters. The free energy differences 
corresponding to the individual panels are (a) AF « —0.87, (b) AF ss —3.81, and (c) AF « —7.28. 


the resulting expression after performing the integration, 
all coefficients depending on A(t/) or A(0) are considered 
as unknown parameters. For example, inserting the po¬ 
tential U(x,X(t)) = k[x — A(f)] 2 /2 from the first line of 
Tab. 1 into Eq. (0, one obtains p(w) = exp[— /3(2w — 

fun 

kAX 2 ) 2 /(8kAX 2 )]/\kAX\ after integration, where Z\A = 
A (tf) — A(0). The constants containing AX are consid¬ 
ered as unknown parameters and marked by A , B and 
C in the table, giving p(w) = Aexp[— (Bw — C) 2 ]. An- 

fun 

other example is the potential U[x , A(t)) = X{t)x 2 /2 from 
the second line of Tab. 1. In this case one obtains from 
Eq. d3]) p{w) = \AXw\- 1 / 2 exp[—/3A(0)u>/ziA]. The result- 

fun 

ing functional form is p(w) = A\w\~ 1 ^ 2 exp[— B\w\], be- 

fun 

cause the exponent —1/2 does not depend on the driving 
AX. 

The first three potentials in Tab. Q] refer to those, 
where exact analytical results for the complete WPD have 
been derived, namely (i) the sliding parabola [T8j, (ii) 
the breathing parabola mm, and (iii) the logarithmic- 
harmonic potential |14j . As a fourth example we added 
the V-potential, where the asymptotic behavior predicted 
by the EN theory could be derived analytically (see Ap¬ 
pendix [B]). For all these potentials, where the asymptotic 
behavior of the WPD is known analytically, the FF con¬ 
jecture on the functional form is correct. 

Let us note that in Eq. m we assumed that AU{x) is 
different from zero, which always is the case for a mono¬ 
tonic protocol (dX(t)/dt strictly positive or negative in 
[0,f/]). Moreover, in Eq. ([31) we not only refer to the 
asymptotic behavior for w —> —oo (relevant for the JE), 
but also to the limit w —»■ +oo. The WPD asymptotic 
behavior in the latter limit is also in agreement with all 
analytically known results. The surprising agreement sug¬ 
gests that the functional form of WPD asymptotic behav¬ 


ior does not depend on the velocity of the driving and 
gives us the motivation to further test the FF conjecture. 

The EN theory is exact in the weak noise limit, i.e., 
up to the first order in thermal energy [see Eq. Q 
in Appendix [X] . The thermal energy determines the free 
diffusion coefficient Dq oc ksT and thus the weak noise 
regime means that the internal dynamics of the system is 
slow. For driven systems, the slowness should be compared 
with the speed of the driving. In this respect, a slowly 
driven system can be obtained by either fixing the driving 
speed and increasing the temperature or vice versa. 

As discussed in the Introduction, we first provide fur¬ 
ther compelling evidence that EN theory is indeed valid 
also for slowly driven systems. We find this necessary, 
because (i) confirmation of this theory by analytical re¬ 
sults has so far been done only for the harmonic poten¬ 
tial, where the propagator for the dynamics is Gaussian, 
and (ii) the EN theory for other potentials has been tested 
by BD simulations, which, however, are restricted to rel¬ 
atively fast driving, where the tails of the WPD can be 
determined with sufficient statistics. For slow driving, the 
essential part of the WPD is always well described by a 
Gaussian distribution m, while the asymptotic behavior 
shows up only when going to extremely large |u>| values. 

The compelling evidence is obtained by utilizing the 
exact analytical results for the logarithmic-harmonic po¬ 
tential U(x,X(t)) = —< 7 log(|a;|) + A(t)x 2 /2 P3j, see third 
line in Tab. |1] because in the presence of the logarithmic 
part (g ^ 0) the propagator for the dynamics is no longer 
Gaussian, and the pre-exponential factor depends sensi¬ 
tively on the strength g of the logarithmic part and the 
inverse temperature j3. For slow driving, the asymptotic 
regime is expected to occur at very large negative work 
values with extremely small p(w) values beyond a possi¬ 
ble physical realization, while the main part of the WPD 
can be well described by a Gaussian. According to the 
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reasons given above, we deliberately want to test the EN 
theory in this situation. 

For this test we first calculate numerically the WPD 
asymptotic behavior using the EN theory (the numerical 
method used for solving the EN equations is described in 
Appendix iDl) and fit the form predicted in Tab.[Q p{w) ~ 
\w\~ v exp(— B\w\), to these data. Then we compare the 
fitted coefficients v and B with the results predicted in 
m- Results for the protocol 

m = ( 4 ) 

l + rt 

with k = 1, r = 1, and t.f = 1 are shown in Fig. [T] for 
three different g values. They strongly suggest that the EN 
theory is correct even for slow driving. The parameters v 
and B : determined from fits in the w interval [—350, w+], 
are shown in the insets as functions of the upper bound 
w+ of the fit interval. As can be seen from these graphs, 
v and B converge to the predicted values with decreasing 
w+. It becomes clear also that the convergence to the exact 
asymptotic functional form occurs at very large negative 
work values, where p(w) attains extremely small values 
of order ICE 300 . Irrespective of this slow convergence to 
the asymptotic functional form, it should be noted that 
the data from the EN theory are very close to the exact 
solution of the complete WPD in the whole w interval 
shown in Fig. [l] 

In the following we now assume the EN theory to be 
valid. The FF conjecture is tested against it and against 
BD simulations, which have been carried out by numerical 
integration of the Langevin equation © with the Heun al¬ 
gorithm [20]. WPD determined from the BD simulations 
have been calculated from typically 10 s simulated trajec¬ 
tories. 

3 FF conjecture versus EN theory for 
equilibrium initial condition 

As mentioned in the Introduction, an analytical solution 
of the differential equations of the EN theory can be ob¬ 
tained for the V-shaped potential, which for any mono¬ 
tonic protocol X(t) gives the WPD asymptotic behavior. 
If dX(t)/dt > 0, p(w) has support only for w > 0, and 
p(w) ~ exp(— Bw) for w —> +oo, while for dX(t)/dt < 0, 
p(w) has support only for w < 0, and p(u>) ~ exp (\B\w) 
for w —> —oo, where B = /3Ao/(A/ — Ao). This functional 
form is predicted by the FF conjecture, where, quite sur¬ 
prisingly, in this case also the parameter B is correctly 
given. We have tested these results against BD simula¬ 
tions (data not shown). 

Further tests of the FF conjecture have been performed 
by us for various potentials, the potentials © and © 
given below, the double-well potential x 4 + X{t)x 2 , var¬ 
ious potentials of the form A(f)|x| n , and the potential 
A(t)log(l + x 2 ). These tests were carried out for differ¬ 
ent protocols, the one given in Eq. ©, the linear proto¬ 
col A(t) = k + rt and the exponential protocol A (t) = 



Fig. 2. Asymptotic behavior of the WPD for the potential 
© as obtained from the EN theory (circles) and BD simu¬ 
lation (squares) for the protocol in Eq. @ with parameters 
k = tf = 1 and r = 10. The three lines are fits of the FF con¬ 
jecture p(w) ~ |w|~ 3 A exp(— B\w\) to the EN data for different 
fitting intervals [—120, w+] with u>+ = —7.0, —4.1, and —1.1. 
The inset shows the fitting parameter B as a function of the 
upper limit w+ of the fit interval. The errors bars given in the 
inset, which are smaller than or comparable to the size of the 
symbols, mark 95% confidence bounds of fitting parameters. 
The free energy difference reads AF « —0.60. 

fcexp(— rt). Our findings suggest that for “hard poten¬ 
tials” , with stronger confinement than the harmonic one, 
the FF conjecture fails in general, and gives, as expected, 
only a good approximate description for sufficiently fast 
driving. 

For “soft potentials”, with weaker confinement than 
the harmonic one, there is evidence that the FF conjec¬ 
ture always predicts correctly the main |w| —> oo asymp¬ 
totics of the WPD (exponential part), and in some cases 
even the pre-exponential factor. This is true for all ex¬ 
amples in Tab. [Q which refer either to soft potentials or 
the harmonic one. It is possible, however, that the pre¬ 
exponential factor is not correctly captured by the FF 
conjecture. 

In the following we exemplify these findings for one 
hard and one soft potential, namely 

U(x, X(t)) = X(t)x 4 , (5) 

and 

U{x,X(t)) = X(t)\x\ + \x\ 3/2 , (6) 

with the protocol A (t) given in Eq. ©. 

Figure [2] shows the results of the EN theory (circles) 
and BD simulations (squares) for the hard potential ©, 
which coincide in the asymptotic regime w -C AF. The 
FF conjecture predicts the asymptotic behavior of WPD 
p(w) ~ |ui| -3 / 4 exp(— B\w\). Fits of this functional form 
to the EN data for three fitting intervals with different 
upper bounds w + are depicted by the lines. As the B 
factor, shown in the inset of Fig. [2] shows no tendency to 
converge to a constant with decreasing upper limit w+ of 
the fit interval [—120, rc+] (see also the non-overlapping of 
the three lines in the main part of the figure), we conclude 
that the FF conjecture fails in this case. 
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duced above (comparison of the confinement with respect 
to the harmonic potential). 

Based on the findings in the previous section, an exten¬ 
sion of the FF conjecture to nonequilibrium initial distri¬ 
butions is considered for soft potentials only. A straight¬ 
forward extension is to replace the p eq ( x) in Eq. ©) by 
p 0 (x), i.e., 


p{w) = 
run 


dx 5 [u> — AU(x)\ po(x), w— >±oo. (7) 


Fig. 3. Asymptotic behavior of the WPD for the potential © 
as obtained from the EN theory (circles) and BD simulation 
(squares) for the protocol in Eq. © with parameters k = r = 
tf — 1. The three overlapping lines are fits of the FF conjecture 
p(w) ~ exp( — Bw — C\w\ 3 ^ 2 ) to the EN data for different fitting 
intervals [—36, W+] with upper bounds w+ = —25.3, —14.5, 
and —3.7. The insets show the fitting parameters B and C as 
functions of w+. The errors bars given in the insets, which are 
smaller than or comparable to the size of the symbols, mark 
95% confidence bounds of fitting parameters. The free energy 
difference reads AF « —0.24. 


Figure [3] shows the results of the EN theory (circles) 
and BD simulations (squares) for the soft potential ©), 
which again coincide in the asymptotic regime w <C AF. 
The FF conjecture predicts the asymptotic behavior p(w) ~ 
exp (—Bw — C |iu| 3 ' 2 ). Fits of this functional form to the 
EN data for three fitting intervals [—36, u> + ] with different 
upper bounds w+ are depicted by the lines. The B and 
C factors from the fits clearly tend to constants, see also 
the overlapping of the three lines in the main part of the 
figure. These lines match well also the simulated data. Ac¬ 
cordingly, the analysis provides evidence that the FF con¬ 
jecture in this case predicts correctly the main |u;| —> oo 
asymptotics of WPD. In the asymptotic formula from the 
FF conjecture, the pre-exponential factor equals one. In 
order to test also this prediction, we have further fitted the 
generalized form p(w) ~ |mj| d exp (—Bw — Cju;| 3 / 2 ) to the 
numerical data. In this case we found that the D factor 
from the corresponding fit shows no tendency to converge 
to zero with increasing upper bound u>+ [see Fig. ©a)]. We 
hence can conclude that the FF conjecture fails to predict 
the pre-exponential factor in this case. 


4 FF conjecture versus EN theory for 
nonequilibrium initial condition 

The generalization of the EN theory to a nonequilibrium 
initial distribution po(x) is presented in the Appendix IQ 
We compared the predicted asymptotic behavior from the 


When testing this ansatz, we indeed found functional forms 
predicted by Eq. 0 to fit the data of the EN theory. How¬ 
ever, this support of Eq. 0 was only obtained, if po{x) 
is “wide” compared to the equilibrium distribution p eq (x) 
in the sense that lim| a .|_ >00 po(x)/p eq (x) = oo. If p 0 (x) is 
“narrow” compared to the equilibrium distribution p eq (x), 
meaning lim| a ,|_ >00 po(x)/p eq (x) < oo, our tests showed 
that Eq. 0 rather than Eq. 0 fits the results of the EN 
theory. From a consideration of protocols with infinitely 
fast driving, this can be understood by studying a one- 
step protocol, where, when starting from the distribution 
Po(x), the potential is suddenly changed by AU (x) at some 
time instant f* > 0 in the interval ]0, tf]. If the propagator 
of the dynamics of the position is known, one can calculate 
the WPD analytically for this one-step protocol by replac¬ 
ing po(x) by the position distribution p(x, i*) at time t* in 
Eq. 0. Investigating potentials with explicit solutions of 
the underlying Smoluchowski equation, we indeed found 
that for wide po(x) the results were in agreement with 
Eq. 0, while for narrow po(x) the results were in agree¬ 
ment with Eq. 0. 

To sum up, the FF conjecture for the non-equilibrium 
initial distribution suggests that for wide initial conditions 
the weights of the trajectories corresponding to large |«>| 
values are determined solely by the initial distribution and 
the potential. On the other hand, for narrow initial condi¬ 
tions the FF conjecture suggests that these weights are de¬ 
termined rather by the evolved distribution at some time 
t* > 0, which in turn yields the same WPD asymptotic 
behavior as the equilibrium distribution. 

To exemplify our findings, we present results here for 
the case of the breathing parabola U (x, A(t)) = A(t)x 2 , us¬ 
ing the protocol A (t) in Eq. 0 with k = r = tf = 1 and 
two nonequilibrium initial distributions po(x) oc exp(—x 4 ) 
(narrow) and po(x) oc x 2 exp(—x 2 ) (wide). As can be seen 
from Fig.© in the asymptotic regime w <C ( w ), where ( w ) 
denotes the mean work done on the system, the results 
of the generalized EN theory (circles) agree well with the 
simulated data (squares) for both initial distributions. For 
the narrow po(x), the FF conjecture predicts the asymp¬ 
totic behavior p(w) ~ |w| -1 / 2 exp(— B\w\), and for the 
wide po{x), it predicts p(w) ~ Ixul 1 / 2 exp(—A?|w;|). Fits of 
these functional forms to the EN data for three fitting in¬ 
tervals with different upper bounds w+ are depicted by 
the lines. The B factors from the fits, shown in the in¬ 
sets, tend to approach constants for decreasing w+. In 


generalized EN theory for various initial distributions againstparticular the data in Fig. ©b) for the wide po(x) have 
extensive BD simulations and found always a good agree- almost reached a plateau value [see the B scale in the 
ment, both for hard and soft potentials in the sense intro- inset in comparison to the one in Fig. ©a)]. Note also 
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Fig. 4. Asymptotic behavior of the WPD for the breath¬ 
ing parabola U(x,X(t)) = X{t)x 2 for a nonequilibrium ini¬ 
tial distribution po{x) and for the protocol in Eq. (j4j) with 
k = r = tf = 1. The circles refer to the data calculated 
from the EN theory and the squares refer to data obtained 
from the BD simulations. In (a) the initial distribution is nar¬ 
row [po(*) oc exp(— x 4 )] and in (b) the initial distribution is 
wide [po(x) oc x 2 exp(— x 2 )}. The three overlapping lines are fits 
of the FF conjectures, p(w) ~ |w| 1 ' 2 exp(—-B|u>|) in (a) and 
p(w) ~ |w| 1 ^ 2 exp( — f?|ui|) in (b), to the EN data for different 
fitting intervals [—100, w+] with the upper bounds w+ = —7.0, 
—4.1, and —1.1. The insets show the fitting parameters B as 
functions of w+. The errors bars given in the insets, which are 
smaller than or comparable to the size of the symbols, mark 
95% confidence bounds of fitting parameters. The mean works 
are (a) (w) = —0.26, (b) ( w } = —0.49. 


that the three fitting lines in the main parts of Fig. 0a) 
and (b) can hardly be distinguished on the scale of the 
graphs, which gives further support that the FF conjec¬ 
ture predicts correctly the main |u;| —> oo asymptotics 
of WPD. Also in this case we have tested the predicted 
pre-exponential factors |w| ±1 - /2 by fitting the generalized 
formula p(w) ~ \w\ c exp(— B\w\) to the numerical data. 
The results suggest that in the case of narrow initial con¬ 
dition the factor C from the corresponding fit shows no 
tendency to converge to —1/2 [see Fig. 0b)]. On the other 
hand, for the wide initial condition the factor C from the 
corresponding fit clearly converges to 1/2 with decreasing 
w+ [see Fig. 0c)]. We hence can conclude that the FF 
conjecture fails to predict the pre-exponential factor for 
the narrow initial condition, while its prediction is precise 
for the wide initial condition. 

To summarize, the main |«>| —> oo asymptotic behavior 
predicted by the generalized EN theory agrees well with 
simulated data, and there is good evidence that the FF 
conjecture for nonequilibrium initial distributions holds 
for soft potentials. The pre-exponential behavior predicted 
by the FF conjecture may be incorrect. 


5 Summary and perspectives 

We have shown that the functional form of the WPD 
asymptotic behavior for overdamped Brownian motion in 



Fig. 5. Test of pre-exponential factors predicted from the FF 
conjecture. In (a) we used the potential © and the equilibrium 
initial condition p eq (*). The panel shows the fitting parameter 
D from the asymptotic form p(w) ~ \w\ D exp (—Bw — C|w| 3 ^ 2 ) 
as a function of the upper bound w+ of the fitting inter¬ 
val [—36, w+]. Panels (b) and (c) correspond to the potential 
U(x,X(t)) = X(t)x 2 and a nonequilibrium initial distribution. 
In (b) the initial distribution is narrow \po{x) oc exp(—a: 4 )] and 
in (c) the initial distribution is wide [po{x) oc x 2 exp(—a; 2 )]. The 
panels show the fitting parameter C from the asymptotic form 
p(w) ~ \w\ c exp(—-B|ic|) as a function of the upper bound w+ 
of the fitting interval [—100, iw+]. The errors bars given in the 
figures, which are smaller than or comparable to the size of the 
symbols, mark 95% confidence bounds of fitting parameters. 


confining potentials can be obtained often from a simple 
conjecture, in which it is assumed that the respective form 
is independent of the driving velocity. This conjecture was 
motivated by the fact that it is valid for all potentials and 
protocols, where analytical solutions for the WPD have 
been obtained in the literature. In fact, these potentials 
are either the harmonic one or they belong to what we 
classified as soft potentials, with weaker confinement than 
the harmonic potential. 

Tests of the FF conjecture have been performed against 
the predictions of the EN theory. With respect to that the¬ 
ory, (i) we gave further evidence of its validity, in particu¬ 
lar in the case of slow driving, by a comparison of results 
with the exact WPD for the logarithmic-harmonic poten¬ 
tial, (ii) we have obtained an analytical solution for the 
V-potential in addition to the formerly known parabolic 
potentials, and (iii) we have extended the theory to the 
case of nonequilibrium initial distributions of the particle 
position. 

The tests of the FF conjecture indeed showed good 
agreement for soft potentials. For hard potentials, with 
stronger confinement than the harmonic potential, our 
tests strongly indicate that the FF conjecture is not valid 
both for equilibrium and nonequilibrium distributions of 
the initial particle position. In the nonequilibrium case for 
soft potentials, the tests were most convincing for wide ini¬ 
tial distributions, which decay slower than the equilibrium 
one. However, also for narrow initial distributions, decay¬ 
ing faster than the equilibrium one, our results indicate 
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the FF conjecture to hold true. All our new findings were 
also checked by BD simulations whenever it was possible. 

It is surprising that for all analytically known cases 
the FF conjecture is valid. Whether this is a mere coinci¬ 
dence or has a deeper reasoning needs to be clarified. To 
this end, theoretical insight should be gained, why the FF 
conjecture fails for hard potentials. 

With respect to applications, the FF conjecture can be 
useful to guess reasonable functional forms for the WPD 
asymptotic behavior, which, as pointed out in the Intro¬ 
duction, are needed, when measured WPD data have to 
be extended to the tail regime for use of the JE. A quite 
general form of the WPD asymptotic behavior was sug¬ 
gested in [5] and, based on it, a “Jarzynski estimator” for 
better determination of free energy differences in experi¬ 
ments was proposed. In cases where it is possible to guess 
functional forms of the driving potential, the FF conjec¬ 
ture can be more specific and utilized to develop improved 
Jarzynski estimators, or to choose a most appropriate one 
among the possible driving potentials. If an appropriate 
potential is found, a modeling of the process by BD sim¬ 
ulations and an application of the EN theory can be used 
to identify the parameters in this potential. 

In our work here we considered overdamped Brownian 
motion in one dimension under monotonic protocols. It 
would be interesting to check in future studies, whether ex¬ 
tensions of the FF conjecture (and also of the EN theory) 
to non-monotonic protocols, to Brownian motion in higher 
dimensions, and to the underdamped case are possible. In 
fact, the FF conjecture has been seen to be valid also for 
a specific case of underdamped diffusion in a breathing 
parabola El- 
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A EN theory 

According to mm, the WPD asymptotic behavior for 
overdamped Brownian motion of a particle in a time- 
dependent confining potential U(x , A (t)) [Eq. Q], for the 


initial particle position in equilibrium with distribution 
Pe q (x) = exp[—/3[/(x,0)]/%q, reads 


p(w) = 


\[2N exp(— j3S) 


^eq 


(det A)(V r, |A” 1 |V r/ ) 


=(i + o(/r 1 )), (8) 


where 


M = exp 



det A = 2 (V t ''xt f +Xt f ), 

{V'\A~ 1 \V') = /W, 

JO 


and 


f tf 

S = Vo + / dt 
Jo 


\{vt + v') 2 + \v 


q 

2 W: 


(9) 

( 10 ) 

( 11 ) 


( 12 ) 


where yt stands for the most probable trajectory which 
yields the work value w (yt = dyt/dt). In all the above 
formulas the mobility p [cf. Eq. (JT])] was set to unity, and 
we have abbreviated the time dependence by the lower 
index. The function V t is defined as V t = U(y t ,X(t)), 
V t = dU(x,X(t))/dt\ x = Vt , Vt = dU(x,X(t))/dx\ x = Vt , and 
analogously for higher derivatives. 

The “optimal trajectory” y t minimizing the action S 
and the auxiliary functions ipt and Xt are obtained from 
the following system of three second-order ordinary differ¬ 
ential equations 


yt = (q- I)v; + VM, (13a) 

i>t = (Vt 2 + V'V'" - (1 - q)Vn^ t - VI , (13b) 

Xt = (V" 2 + VX" - (1 - q)V t ")xt, (13c) 


with the boundary conditions 


2/o = Vq , y tf = -Vl f , 

ipo = Vo'ipo , ip tf = —V" f ip tf , 

Xo = 1, Xo = Vq . 


The constraint 


w= dtVt 

Jo 


fixes the q value in Eqs. (fl2l) . (I13bl) . and « m - 


(14a) 

(14b) 

(14c) 

(15) 


B EN theory applied to V-potential 

For the potential U(x,X(t)) = A(i) |ac| the right hand side 
of Eq. © can be calculated analytically. As we show be¬ 
low, for large enough values of |u)| the optimal trajecto¬ 
ries yt never crosses the origin. Since we are interested 
in the tail behavior |w| —> oo, this means that the non- 
analytical nature of |x| at the origin does not influence 
the asymptotic behavior of the WPD and in the Eqs. @- 
m, one can safely assume that V" = V"' = 0 etc. for 
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all x. Then the differential equations for the functions ipt 
and Xt, which determine the pre-exponential factor in (|8jl 
become independent of yt and thus of w also. 

Due to the symmetry of the problem, for each work 
value w two optimal trajectories yt symmetric relative to 
the origin exist. In the following, we will solve the equa¬ 
tions presented in the preceding section assuming that 
yt > 0. The trajectories yt < 0 can be incorporated sim¬ 
ply by including the factor 2 in the prefactor of the WPD. 
The optimal trajectory y t follows from the equation 

m = {q- i)A t , 

which can be easily solved. The two unknown integration 
constants together with the auxiliary variable q are deter¬ 
mined from the formulas 


f tf 

Vo = X 0 , 2 )t f = -A t f , / dtV t = w . 

Jo 

The solution is (for monotonic driving, where A t f yf Aq) 


Ut — ft + 


Xt f — Aq 


w, q = - 


2An 


A t f — Aq 


(16) 


where ft is a function independent of w. Due to the work 
definition m, the term w/(Xt f — Xo) is for a monotonic 
driving always positive. The function f t is finite and thus 
there always exists some \w\ large enough that the whole 
function y t is positive for an arbitrary t. The action (1121) 
can be rewritten as [12] 


S ~ + 2 (^/ + Vo) ~ j (2 h f Vt f +2/0^0) 

+ i J* f dt Vf{v; - y t V”) + J* f dt y t v; . (17) 

After inserting the optimal trajectory (1161) into this for¬ 
mula one obtains 


S — /r + 7(At/2/t/ + Aq?/o) + 


dtXi 


1 -q 


= h 


tf 


Xt f ~ Aq 


w, 


where h t is a another function independent of w. Accord¬ 
ingly, the WPD fulfills the asymptotic relation 


p(w) 


exp xf f X -\ 0 w ) > w +°° > At > 0 , 
exp Ao ^ w'j , w — > —00 , At < 0 . 


(18) 


C EN theory for nonequilibrated initial 
condition 

The EN theory m can be generalized to a nonequilibrium 
initial condition po(x) = exp[— /3r(x)\/Zo of the particle 
position. The derivation is straightforward and we present 
only the results here. 


The asymptotic form (0 remains valid, but the action 
is modified, 


rt f 

s = r+ dt 

Jo 


i fa 


V 


/\2 


\ v 


q 

~2 W ’ 


(19) 


where T = r(yo). Also the auxiliary functions yt, ipt 
and \t for determining S from Eq. m and detA and 
<W'|^4- 1 |V' , > fromEqs. IfTTTl) and (fill) , respectively, become 
modified. Their evolution equations m remain the same, 
but the boundary conditions now are 


yo = 2 r' - v', 

II 

(20a) 

tp 0 = (2 r" - v 0 >o , 

II 

--V 

(20b) 

Xo = 1 , 

xo = 2 r" - v 0 ". 

(20c) 

D Numerical procedure 

used for solving 

EN 


equations 

In order to solve the EN equations (Appendices [A] and 
ED, we have adopted the numerical procedure “bvp4c” 
implemented in MATLAB and previously used by Nick- 
elsen da- Bvp4c is a finite difference code, which solves 
two-point boundary value problems for ordinary differen¬ 
tial equations with possible further unknown parameters. 
It implements the three-stage Lobatto Ilia formula EH- 
This is a collocation formula and the collocation polyno¬ 
mial provides a continuous solution with continuous first 
derivative that is fourth order accurate uniformly in [a, b\. 
Mesh selection and error control are based on the residual 
of the continuous solution. 

In the code, an initial guess of the solution is used to 
iteratively find the correct solution. In our implementation 
we guessed the solution for some small absolute value of w, 
and use the corresponding solution as the initial guess for 
a new, slightly changed, work value. The corresponding 
solution is then again used as the guess for another work 
value and so on until a given final value of w is reached. 
Problems arise when the solutions of EN equations for 
two close work values differ significantly. Then the numer¬ 
ical procedure collapses, because it is not able to converge 
from a distant initial guess to the solution. In our analysis 
we have encountered such numerical problems when using 
very soft potentials such as U(x, X(t)) = X(t) ln(l + x 2 ). 
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